Quantum discord dynamics of two qubits in the single-mode cavities 
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The dynamics of the quantum discord for two identical qubits in both two independent single- 
mode cavities and a common single- mode cavity are discussed. For the initial Bell state with 
correlated spins, while the entanglement sudden death can occur, the quantum discord vanishes 
only at discrete moments in the independent cavities and never vanishes in the common cavity. 
Interestingly, quantum discord and entanglement show opposite behaviors in the common cavity, 
unlike in the independent cavities. For the initial Bell state with anti-correlated spins, quantum 
discord and entanglement behave in the same way for both independent cavities and a common 
cavity. It is found that the detunnings always stabilize the quantum discord. 

PACS numbers: 03.65.Ud, 75.10.Jm, 03.67.Mn 



I. INTRODUCTION 

Quantum entanglement, originated from nonlocal 
quantum correlation, is fundamental in quantum physics 
both for understanding the nonlocality of quantum me- 
chanics [T] and plays an important role in almost all ef- 
ficient protocols for quantum computations and commu- 
nications [2]. Due to the unavoidable interaction with 
the environment, an initially entangled two-qubit system 
becomes totally disentangled after evolving for a finite 
time. This phenomena is called entanglement sudden 
death (ESD) [3 1 and has been recently demonstrated ex- 
perimentally (4|. However, the entanglement may fail to 
capture the existence of the quantum correlation in some 
mixed separate states, in which the entanglement is con- 
sidered not a good measure 0, Q . 

Recently, a new kind of the quantum correlation, quan- 
tum discord (QD) has attracted a lot of attentions j7|]. It 
provides the alternative route for measurement, which 
is present even under separable states Q. The defini- 
tion of the QD can be interpreted as the difference of 
the total quantum information of the two sub-systems A 
and B before and after the local operation on the one of 
them. The QD has been proved as a good measure of the 
non-classical correlations beyond entanglement. Further- 
more, the QD has been indicated as the source to speed 
up the quantum computations ^8, 9J. 

Some works have been devoted to the QD dynamics 
of two qubits coupled to Markovian jlo|, and non- 
Markovian [T2j environments. The comparisons with en- 
tanglement dynamics have been also performed. How- 
ever, the relevant study on two two-level atoms (qubits) 
coupled to independent or common single-mode cavities 
without dissipations has not been found in the litera- 
ture, to the best of our knowledge. We believe that the 
QD dynamics in these qubit systems is also of fundamen- 



tal interest. In addition, some essential pictures can be 
clearly described and unfolded in the framework of the 
simple model where the exact solutions are available. Ac- 
tually, the entanglement dynamics for two independent 
Jaynes-Cummings (JC) atoms has been well studied pre- 
viously [IMl. The ESD was observed obviously from 
the initial Bell states with correlated spins. This feature 
would prevent the application of the entanglement as ba- 
sic resource for quantum information processing. What 
is the consequence for the QD in this kind of the qubit 
system?. It is just the main topic of the present study. 

In the present paper, we will study the QD dynamics 
for two identical qubits in both two independent iden- 
tical single-mode cavities and one common single-mode 
cavity. Comparisons with the corresponding pairwise en- 
tanglement, i.e. concurrence, are also given. The paper 
is organized as follows. In Sec. II and HI, we derive the 
time dependent QD in these two systems if initiated from 
two typical Bell states. In Sec. IV, the results are given 
and discussions are made. The conclusion is presented in 
the last section. 



II. QD IN TWO IDENTICAL 
JAYNES-CUMMINGS ATOMS 

The Hamiltonian of two identical Jaynes-Cummings 
atoms is shown as 
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, y, z) is the Pauli operator of the atom 
A(B), shown as (7, = I t > ( t I - I i ) ( i I and a, = I t ) ( i 
I -I- I |)( t I, with I t )(| I >) the excited (ground) state of 
the two- level atom, a\b^) and a{b) are the creator and 
annihilator of the cavity A (B), respectively, A and u 
are the frequencies of the atom and the cavity, g is the 
atom-cavity coupling strength. Here we set Ti = 1 and 
the detunning 5 — IS. ~ lo. 
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We first study the evolution of the QD initiated from 
the Bell state with anti-correlated spins, which has the 
following form 



*Bcii> = sina|;t) +cosa|t;). 



(2) 



Initially, the vacuum state of the cavity is considered, so 
the initial state of the whole system can be written as 



|*(0)) = (sina|;t) +cosa|t;>)®|00) 
^ sina|;tOO> +cosa|UOO). 



(3) 



The time dependent wave function can be generally ex- 
pressed as [1^] 



where the coefficients are 
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-|-Be"'^-*)cosa 
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X2 


= {Ae 




+ Be" '^-*) sin a 




X3 


= C(e 


— iA-i- 1 


-e^'^-*)cosa 




X4 


= C(e 


— iX-\- 1 


- e" '^-*) sin a. 





with the eigenfrequencies as 



2^ : 



G2 



(6) 



here G = 2g. The auxiliary parameters are shown as 

1 S 
2 
1 

2^' 
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A 
B 
C 



V(52 + G2 
5 



(7) 



2V<52 + G2 



The pairwise density matrix from Eq. under the 
standard basis |it)i lti)i Itt)} is thus expressed 

by tracing the freedoms of the cavities pAB{t) = 
Tr,,4p(t)} = Tr,,„{|*(t))(*(t)|}, 



PAB{t) = - 



( |a;3p + |a;4|2 

|x2p x\x-i. 

X2X1 \x\ 

V 00 



2 

, 



(8) 



With this density matrix, the routine to derive the QD 
is formally given in the Appendix A. The von Neumann 
entropy for two atoms in Eq. (jA2l) is given by 

(9) 



S{A,B) = -(|x3p + |x4niog(|a:3p + |a;4n 
2 + |x2p)log(|a:ip + |a;2n, 



and the sub-system entropies in Eq. 

S{A) = -(l-lxip)log(l 
~\xx\^\ogW\\ 

S{B) = -(l-|x2|2)log(l 
-|a;2plog|x2p. 



3]) are shown as 
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From the Appendix A, one can find that the expres- 
sions of the elements in Eq. (|Aip are 



t;+ = |x3p + |a;4p, w_ = 0, w= \x2^ , 



y = xxx*^, u = 0. 



(12) 



Moreover, 
Xi 



Xl.+ = (|x3p + |x4nC0s2( 



Ixif COS^ ( 



= '^\x,\'\x,\'. 



|a;2psin^( 



and 



a;3p + sin^ I 



X2,+ - 

X2.- = \xi\''sm'9, 

,2 a 



|a;2p cos^ ( 



\Y2\^ 



sm 



\X1\'\^2\'. 



Therefore these parameters are independent of </>. It fol- 
lows that we can search the minimum of the conditional 
von Neumann entropy by only varying in the regime 
[0,7r/2]. Following the procedures outlined in Appendix 
A, we can finally derive the quantum discord. Since a is 
limited to (0,7r/2) , it can be numerically checked that 
9 — tt/A corresponds to the minimum of the conditional 
entropy in the following calculations. The minimum of 
the conditional von Neumann entropy reads 



S{A\Il^ 



Ve log Ve 



(13) 



where 

77± = {1±[(1 - 2\x,\Y + Mmn'/'}/2. (14) 
As a result, the quantum discord is finally given by 

V = -(1- 1:^21') log(l-|a:2n-|x2p log |X2P 
+ {\x3\^ + \x^\^)\og{\xs\' + \xi\^) 

+ (|xip-Mx2ni0g(|xiP + |x2n 

-^?7,log??,. (15) 

For later use, we also list the expression for concurrence 
derived in Ref. [ISj as 



GAB{t) = I sin 2a I [1 -4(7^ sm^{y/5^ + GH/2)]. (16) 

Next, we consider the Bell state with correlated spin 
as the initial atomic state, which is 



l^-^j'^^i) -sinalU) +cosa|n). 



(17) 



Including the initial vacuum cavities, the wave function 
of the whole system can be expressed as 

\^){t) - xi|n00)+a-2|Ull)+a^3|n01) 

+a;4|;tlO) +X5IUOO), (18) 



3 



where the coefRcients are 



Xi 










(19) 


X2 












X3 






+ Be-'^- 




cos a, 


X4 


= C(e-'^+*- 




+ Be~'^- 




cos a, 


X5 


= sin a, 











The eigenfrequencies and the auxihary parameters are 
the same as those in Eqs. (HJ and ([7]). Then under the 
standard basis \it): Iti): ltt)}i the pairwise density 

matrix is shown as 



/\X2\' + \X5\' 



PAB{t) = 





\X4\' 



V 



\X3 








xix^ 

Hence, the joint von Neumann entropy is derive as 

S{A,B) = -|x3plog|a:3p-|a:4plog|x4p (21) 



£=± 



fie log fie 



where 



{{\xi\' 



\X2\' + \X5?) 



(22) 



±Vi\^2\' + ksP - \xi\^r + 4|xiP|x5P}/2 
And the sub-system entropy can be derive as 

S{A) - -{\X2\^ + \X,\^ + \x5\')\0gi\x2\' + \X4\' 

-{\x,\^ + \xs\^)logi\x,\' + \x3\'), 
S{B) = ~i\x2\' + \xs\^ + \x5\')\ogi\x2\' + \x3\' 

-(|xi|2 + |a;4|2)log(|xip + |x4n. 



x,\') 
(23) 

(24) 



Similar to the above Bell state with anti-correlated spins, 
if we focus on aG(0, 7r/2), 6 = 7r/4 also corresponds to 
the minimum of the conditional von Neumann entropy. 
Hence, the minimum of the conditional entropy is given 
by 



S{A\n^) = -^rjJogv,, 



(25) 



where 



V± = {1±V(1 - 2|a;i|2 - 2|X3|2)2 +4|xi|2|a;5|2}/2.(26) 

As a result, the quantum discord can be derived from 
Eqs. (1231), (IMl), and (HI]) as 

V = S{B)~S{A,B) + S{A\n^). (27) 

The concurrence in this case has been also derived pre- 
viously [31, and is also collected here 



CAB(i) = max{0,/(0} 



fit) = [l-iC^sm^{VS^ + GH/2)] 



(28) 

[I sin cos^ a] . 

Specially at resonance (6 = 0), the entanglement sudden 
transition occurs only for a < ac, where ac = it /A. 



III. QD IN TWO IDENTICAL QUBITS IN ONE 
COMMON SINGLE-MODE CAVITY 



The Hamiltonian of two identical qubits interacting 
with one common single-mode cavity reads 

k=A,B 

where and a are the creator and annihilator of the 
common cavity. Actually, it is just the N = 2 Dicke 
model [l^. The detunning is also set as (5 = A — 

If the initial atom state is selected as the Bell state with 
anti-correlated spins, we can obtain the time dependent 
wavefunction as 



i*(i)) - xiiw) + X2\m + xsiui); 



(29) 



with 



Xi 


= (no|e- 


-iHuN2t 




X2 


= (;to|e- 


-iHr)N2t 




X3 




-iHuN2t 





Then the pairwise density matrix under standard basis 
{|U),|;t),|n),in)} is given by 



PAB[t) 



|x3p 

|X2P xlX2 

XiX*2 

0, 



(30) 



For the resonant case (5 = 0), we specify the coefficients 
of the wavefunction in Eq. as 



Xl 
X2 
X3 



COS a 
2 

cos a 



(cosV2Ai + 1) 
(cosV2Ai- 1) 



sma 
2 

sin a 



(cosV2At- 1) 
(cos\/2At-|- 1) 



V2 



(cos a + sin a) sin \/2Xt. 



After the numerical checks, we find that 6 = w/A cor- 
responds to the minimum of the conditional entropy at 
arbitrary time. Hence the QD is described as 

V = -(l-|x2niog(l-|a;2n-|a;2plog|x2p 

+ (1X3^ l0g(|x3|^) + (1 - 1X3^ l0g(l - |X3|2) 
y^77fl0g77f, (31) 



^ 77e 10g77e, 



with 



ry± = {1±[(1 - 2|xi|2)2 + A\xiX2\^]'^'}/2. (32) 



Besides, the concurrence of the two atoms can also be 
given as 



CAsit) = 2max{G, |xiX2|}. 



(33) 
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FIG. 1: (Color online) Resonant dynamics of QD and concur- 
rence of two identical JC atoms with the initial atomic Bell 
states with anti-correlated spins (a) and correlated spins (b). 
ai = 1.0. The blue solid line and the green dashed line are 
corresponding to QD and concurrence respectively. 



If starting from the initial Bell state with correlated 
spins, the wavefunction can be expressed as 

= xiino> + a;2i;;2) + a;3ini> + a;4i;ti> 
+x5i;;o). (34) 

Hence the pairwise density matrix can be derived as 



/|a:2p + |x5|2 



PAB{t) 








1 2:4 1 ^3X4 
x:ixl |a;3p 



(35) 











The coefficients at resonant condition are shown as 
cos a 



Xl 

X2 

X3 

X4 
X5 



.e-*"*(2 + cos VQXt), 



V2 



cos a 



3 

-i cos a 



(cos V&Xt — 1), 



V6 

—i cos a 

sin a. 



e-"^* sinVeXt, 



FIG. 2: (Color online) Off-resonant dynamics of QD for two 
identical JC atoms with the initial atomic Bell states with 
anti-correlated spins (a) and correlated spins (b) for different 
detunnings S — g,2g, and 5g. ui = 1.0. 



While for the concurrence at resonance, we know that 
tlie pairwise density matrix in Eq. p5|) has the form as 



PAsit) 



/ v+ u* 

y y 

y y 

\u V- 



Then we can derive it as 



CAsit) = 2max{0, |a;iX5| - |a;3p, 

\x3?-\xiW\x2\' + \x5n- 



(36) 



(37) 



Then we can obtain the QD numerically. 



From the definition, we find that there exists a critical 
bound for a. The ESD happens only for a < ac- The ac 
is determined by 

tan^ Qfc — 4 tan + I — 0, (38) 

resulting in ac = arctan(2 — -\/3) — tt/12. 



IV. RESULTS AND DISCUSSIONS 

First, we compare the QD with the concurrence in the 
two identical JC atoms with two initial atomic states, i. e. 
the Bell states with anti-correlated spins and correlated 
spins, for zero detunnings. The results are collected in 
Fig. 1. The evolution of both QD and concurrence for 
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Gt/:it 




FIG. 3: (Color online) Resonant dynamics of QD and con- 
currence of two qubits coupled to a common cavity with the 
initial atomic Bell states with anti-correlated spins (a) and 
correlated spins (b). uj = 1.0. The blue solid line and the 
green dashed line are corresponding to QD and concurrence 
respectively. 



the initial Bell state with anti-correlated spins display 
similar behavior, as seen in Fig. 1(a). Yonac et al. 
has shown that the ESD only occurs in the initial atomic 
Bell states with correlated spin, where the entanglement 
can fall abruptly to zero and vanish for a period of time 
before revival. It is interesting to note from Fig. 1(b) 
that during the period of ESD, QD becomes small but 
is always finite, except vanish at discrete moments t = 
(2fc-|-l)7r/G, (fc = 0, 1,2,...). 

Then, we show the effects of the detunnings on the QD 
in independent cavities in Fig. 2, starting from these two 
Bell states. Interestingly, the amplitude of oscillation of 
the QD as a function of time is suppressed monotonically 
by the detunnings for both initial Bell stats. More im- 
portantly, the zeros of the QD at discrete instants shown 
in Fig. 1(b) disappear with the finite detunnings. 

Next, we compare the QD with the concurrence in the 
two identical qubits coupled to the common cavity with 
two initial atomic Bell states. The results for zero de- 
tunning are presented in Fig. 3. For the initial atomic 
Bell state with anti-correlated spins, similar behaviors for 
both QD and concurrence are observed. For the initial 
atomic Bell state with correlated spins, one can find from 
Fig. 3(b) that the ESD can occur, but the QD never van- 
ishes. Interestingly, QD and entanglement show opposite 



FIG. 4: (Color online) Off-resonant dynamics of QD for two 
qubits coupled to a common cavity with the initial atomic 
Bell states with anti-correlated spins (a) and correlated spins 
(b) for different detunnings 5 = g, 2g, and 5g. u) — 1.0. 



behaviors. Especially, during the period of ESD, the QD 
always becomes larger, in sharp contrast with that ob- 
served in the independent cavities (c.f. Fig. 1(b)). 

Note that the critical parameter Uc = tt/12 below 
which the ESD can occurs in the common cavity is 
smaller than ac — vr/4 in independent cavities. The 
instant vanish of QD is absent in the common cavity, 
implying that the common cavity enhance the QD. So, 
it is suggested that the quantum correlation in the com- 
mon cavity is stronger than that in independent cavities 
in some sense. 

The effect of the detunnings on the QD of two qubits 
coupled to a common cavity is also studied. As shown in 
Fig. 4(a), for the initial Bell state with anti-correlated 
spins, the amplitude of oscillation of the QD as a function 
of time is also suppressed monotonically by the detun- 
nings and little bit larger than that in two cavities (c.f. 
Fig. 2(a)). While for the initial atomic Bell state with 
correlated spins with a = 7r/24 where ESD can occur, the 
oscillation of the QD for two qubits in the common cavity 
is suppressed considerably with detunnings, as shown in 
Fig. 4(b). Especially, for large detunnings S = 5g, the 
QD remains almost unchanged. In this case, we find that 
the components of the Bell stats show slightly variation 
with time for large detunning, due to the fact that large 
detunnings prevent the hopping for photons between dif- 
ferent atomic levels to certain degree. 
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V. CONCLUSIONS 

In this paper, the QD dynamics of two qubits in both 
independent and common cavities are investigated. The 
comparisons with the entanglement evolution are also 
performed. For the initial atomic Bell state with anti- 
correlated spins, the QD and entanglement show the sim- 
ilar behaviors for both cavities. But for the initial atomic 
Bell state with correlated spins, the QD and entangle- 
ment behave in a remarkably different way. The ESD 
may occur for both cavities, but the QD never vanishes 
suddenly. For the independent cavities, the QD vanishes 
only at discrete instants and can be lifted with finite de- 
tunnings. In the common cavity the QD is always finite. 
Especially, the QD and entanglement display an oppo- 
site behavior in the common cavity, different from those 
in independent cavities. The detunnings play important 
role on the QD dynamics. It always stabilizes the QD, 
which could be helpful in the real applications of the QD 
as the better resource in quantum information science 
and quantum computing. 
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Appendix A: Derivation of the quantum correlation 

In the present paper, the general pairwise density ma- 
trix under the standard basis {||t) , lit) , Itl) Jtt}, is 
shown as 



PAsit) 



f v+ u* 

w y* Q 

y X 

\u w_ 



The von Neumann entropy of the two atoms is 

S{A,B)^ J2 i^k,e\0gnk,e, 



(Al) 



(A2) 



e=±,A:=l,2 



with 



{v+ + v^)±y/{v+ - + 4|u|^ 



n 



2.± 



{w + a;)± {w — x)'^ + A\y\- 



The reduced sub-system density matrices for A and B 
are obtained as 

PAit) = iv++w)\i)^(i\ + ix + v_)\t)^{tl 

PB (t) =_{v++x)\l)s{l\ + {w + V-)\t)B{tl 

Hence, we derive the corresponding von Neumann en- 
tropies as 

S{A) = -{v+ +w)\og{v+ +w) 

-{x + v^) \og{x + u_ ) , (A3) 

S{B) = -{v+ + x)\og{v+ + x) 
— {w + V-) \og{w + 



While for the conditional density matrix Paih^, projec- 
tion basis are considered as 



|$i)^ - cos0\i)^ + e''^sme\t)s, 
\^2)b = e-'^sme\i)g-cose\t)B. 

The conditional density operator is expressed as 

PAln^ = ^kPAB^k /Pk, 



where Hf = lA®\'^k) B{'^k\ and pk TrAs{pA|nf}- 



Specifically under the projections in Eq. (jA4l) . 

Pa\U- = \'^k)B{'^k\®{\i)A{i\XK+ + \^)A{^\Xk,- 

+li)A(tin + it)A(tin*}/pfe- 

For fc = 1, we show 

= v+ cos^ 6 + w siv? 9, 

ATi,- = X cos^ +V- sm^ 9, 

Yi = (?;*e~^'^ -f M*e*'*)sin6lcos6'. 

For fc = 2, 

X2.+ = v+ sin^ 9 + w cos^ 9, 
X2^- = X sin^ 9 + V- cos^ 9, 

Y2 = -(y*e"''^-f w*e^'^)sin6'cos6i. 

Then the eigenvalues of the conditional density matrix 
reads 

Vk,± - ■^{{Xk.+ + Xk,-)±[{Xk,+ - Xk,-f 
2pk 

+4|rfc|2]i/2}. (A4) 
The conditional von Neumann entropy is described as 
S{A\Il^) = ^ -pfcTr^{p^|nBlogp^|n-} (A5) 
= - ^Yl PkVk,e{9,(j)) log r]k,ei9,(l)). 

e=± k=l,2 

As a result, the quantum discord can be obtained by 

•D = min{S'(B) - S{A, B) + S{A\U^)}, 
and the classical correlation is given as 

C = max{5(A) - 5(A|n^)}. 
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